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The phenomenon of matter wave interference lies at the heart of quantum physics. It has 
been observed in various contexts in the limit of non-interacting particles as a single parti- 
cle effect. Here we observe and control matter wave interference whose evolution is driven 
by interparticle interactions. In a multi-path matter wave interferometer, the macroscopic 
many-body wave function of an interacting atomic Bose-Einstein condensate develops a reg- 
ular interference pattern, allowing us to detect and directly visualize the effect of interaction- 
induced phase shifts. We demonstrate control over the phase evolution by inhibiting interaction 
induced dephasing and by refocusing a dephased macroscopic matter wave in a spin-echo 
type experiment. Our results show that interactions in a many-body system lead to a surpris- 
ingly coherent evolution, possibly enabling narrow-band and high-brightness matter wave 
interferometers based on atom lasers. 
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Matter wave interference has been observed as a single particle effect for electrons 1 , neutrons 2 , 
atoms and molecules 3 . Macroscopic matter wave interference was first directly observed in the case 
of two independent atomic Bose-Einstein condensates (BEC) that were brought to overlap 4 . This 
experiment validated the notion of the BEC as a macroscopic matter wave and coined the expres- 
sion of the atom laser in analogy to the laser for the case of photons. Matter wave interferometers 57 , 
in particular for applications to precision measurements, are typically operated in the dilute sin- 
gle particle limit 810 to avoid particle-particle interactions. Atom interferometers based on Bose- 
Einstein condensates (BEC) are expected to benefit from the extremely low momentum spread, 
the exceptional brightness, and the low spatial extent of the BEC 11 , but they readily enter the non- 
linear matter wave regime as a result of the interaction-induced mean field potential. A possible 
solution is to operate BEC-based interferometers in the non-interacting limit 12 - 13 by exploiting the 
cancellation of the scattering phase shift near a scattering resonance. This condition, however, 
is difficult or impossible to fulfill for most atomic species. In the present work we demonstrate 
a BEC-based multipath atom interferometer where the dynamics is dominated by interaction- 
induced phase shifts. We realize the multipath interferometer by loading an interacting BEC into 
an optical lattice potential along one dimension, coherently splitting the BEC into several parts that 
are then each subject to different linear and nonlinear phase shifts. The linear phase shifts due to 
the gravitational force lead to the the well-known phenomenon of Bloch oscillations 1415 , whereas 
the interaction-induced nonlinear phase shifts cause the macroscopic wave function to first spread 
in momentum space as a function of time and then, surprisingly, to exhibit high-contrast inter- 
ference. We demonstrate a high degree of coherence by reversing the nonlinear phase evolution, 
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thereby refocusing the BEC momentum wave function. By application of an external potential 
we cancel the dominant mean-field contribution to the phase evolution and become sensitive to 
beyond-mean-field effects. A crucial ingredient of our experiments is the capability to tune a, the 
atomic scattering length which determines the strength of the interaction, by means of a Feshbach 
resonance 16 . In particular, a can be switched to zero to stop the interaction driven part of the evo- 
lution in the interferometer or to perform high resolution wave function imaging in momentum 
space. 

1 Phase evolution in the matter wave interferometer 

Our interferometer consists of a BEC in adiabatically loaded into a ID optical lattice potential 
with a superimposed harmonic trap, as illustrated in Fig. la. In the tight-binding regime, it is 
convenient to write the macroscopic wave function of the condensate, ^, in a basis 17 of wave 
functions ^j{z, r±) centered at the position Zj = jd of the individual lattice sites j, ^{z, r±,t) = 
Y^j Cj(t)^j(z,r±). Here, z is the coordinate along the (vertical) lattice direction, r± is the trans- 
verse coordinate, d is the distant between adjacent lattice sites and Cj(t) are time-dependent com- 
plex amplitudes. 

After the BEC is loaded into the lattice, we tilt the lattice potential by applying a strong 
force F along the lattice direction. In the limit Fd 3> J, where J is the tunneling matrix element, 
tunneling between lattice sites is inhibited. The on-site occupation numbers \cj\ 2 are then fixed, 
and we can write cj(t) = Cj(0)e l( ^ j ^\ where the phase (f>j(t) evolves in time according to the local 
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potential at each specific lattice site 21 , 

h 9 ^ = Fdj + vf ap + ^ c 

= Fdj+p tT j 2 -a- mtJ 2 . (1) 

Here, the total potential at each lattice site j consists of three terms. The applied force leads to a 
term linear in j and causes Bloch oscillations 14,15 with angular frequency Fd/h. The second term 
comes from an optional harmonic confinement, where f3 tr = muo 2 r d 2 /2 characterizes the strength 
of the confining potential and uo tT is the trap frequency. Atom-atom interactions give rise to a 
third term, the local chemical potential fi l ° c , which depends on the scattering length a and the site 
occupation number as 17 fi l ° c oc ^Ja\cj\ 2 . When the BEC is loaded in the Thomas-Fermi regime, 
as is done here, its initial value can be calculated in a simple way. The density distribution will 
be such that the local chemical potential mirrors the trapping potential used when loading into the 
lattice, n\ oc = /x — \/ trap , with /j, being the (global) chemical potential of the BEC. We then initially 
have n l ° c = a int j 2 , where a int = muo 2 d 2 /2 and uo\ Q is the trap frequency during loading. 

The phase terms proportional to j 2 lead to a nonlinear phase evolution and to a dephasing of 
the lattice site. This results in a time- varying interference pattern of the macroscopic matter wave, 
as we will demonstrate below. The key in our experiments is that we have full control over these 
nonlinear terms, not only over (3 tI via the external trapping potential, but also over the interaction 
term characterized by a int , both via the initial density distribution, and, more importantly, via the 
scattering length a. By tuning the scattering length 16 from its initial value a to a', we can ramp 
a int to a new value a[ nt , which can in particular be set to zero for a = 0. Nonlinear phase terms 
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for matter waves are well known in single particle quantum mechanics. They play an important 
role for matter wave Talbot interferences 5 ' 18 and can be visualized in terms of so-called matter 
wave quantum carpets 19 . In these contexts, the phase terms arise from propagation. In our case, 
the nonlinear phase terms for a- mt ^ arise from interactions and thus lead to a density dependent 
many-body effect in the multipath atom interferometer. 

In the preceding discussion, we have assumed that the minimum of the trapping potential 
is centered directly over one of the lattice minima. If this is not the case, the trapping potential 
term in equation (1) has to be modified to f3 tr (j — 5) 2 = (3 tr j 2 — 2f3 tr 5j + const., where S G [0, 1] 
describes the offset of the trap center in the ^-direction with respect to the lattice minima, and an 
analogous modification has to be done to the interaction term. This adds a small term linear in j 
and therefore leads to a slight modification of the Bloch oscillation frequency. In our experiments, 
5 is the only parameter that we do not fully control. It is constant on the timescale of a single 
experimental run, but it drifts over the course of minutes as the beam pointing of the horizontally 
propagating laser beam generating the trapping potential is not actively stabilized. 

2 Interaction-induced matter wave interference 

The starting point for our experiments is a BEC trapped in a crossed optical dipole trap and adiabat- 
ically loaded into an optical lattice, as illustrated in Fig. la. The gravitational force acting on the 
BEC is initially compensated using magnetic levitation 16 . We effectively start the multipath atom 
interferometer and hence the evolution of the interacting macroscopic wave function by turning off 
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magnetic levitation and ramping down the vertical confinement created by laser beam L 2 within 
0.3 ms, inducing Bloch oscillations in the lowest band of the lattice. With Fd/h m 2n x 1740 Hz 
and J/h ps 2ir x 40 Hz the on-site occupation numbers \cj\ 2 are fixed to their initial values. After 
an evolution time r, we close the interferometer by ramping down the lattice in 1 ms and directly 
image the (vertical) quasi-momentum distribution in the first Brillouin zone (BZ). The ramp is 
adiabatic with respect to the bandgap and maps quasi-momentum onto real momentum 20 , which is 
measured by taking an absorption image after a period of free expansion. Fig. lb shows absorption 
images of the first Bloch oscillation 14 . The Bloch period is about 0.58 ms and the peaks have a 
root mean square (rms) width of 0.2hk, where k = n/d is the lattice wave vector, thus being well 
separated. 

We study the evolution of the wave function at high resolution in momentum space by taking 
snapshots after extended time-of-flight. As illustrated in Fig. 2a, the BEC wave function spreads 
out in the BZ in about N = 18 Bloch cycles. Then, surprisingly, an interference pattern gradually 
develops at the edge of the BZ and later also becomes visible at the center of the BZ, while the 
number of interference maxima and minima changes as time progresses. Images are taken after 
an integer number of Bloch cycles for cycle phase = 0, corresponding to the first image in Fig. 
lb. The data is acquired with an interacting BEC with the scattering length set to 190 a , where 
a is the Bohr radius, at an initial peak density of n = 4 x 10 13 atoms/cm 3 , occupying about 35 
lattice sites after loading. We can follow the evolution of the interference pattern for more than 
N = 100 Bloch cycles, corresponding to times beyond 60 ms. This is about a factor 10 longer than 
the timescale for the initial broadening. We find that the number of maxima and minima and their 
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location in the interference pattern as measured after fixed evolution time r depend on the initial 
atomic density, on the strength of the interaction, and on the number of occupied lattice sites. We 
also find that the measured quasi-momentum distribution for a given r is reproducible from one 
experimental realization to the next, except that the pattern appears slightly shifted within the BZ 
after several experimental realizations. We attribute this to a drift of S, the offset of the lattice 
minima from the dipole trap center, which leads to a small change of the Bloch frequency as noted 
before. We do not actively stabilize the vertical position of L 2 with respect to the lattice, and hence 
temperature variations in the laboratory slowly change 5. 

We combine two techniques to achieve a high resolution in momentum space and to visualize 
the interference pattern. First, we minimize broadening of the distribution as a result of interactions 
by setting a to zero during the release from the lattice and the subsequent free expansion 12 . In 
addition, we use long expansion times, employing magnetic levitation to prevent the BEC being 
accelerated by gravity and falling out of the field of view. Fig. 3 shows how the contrast emerges 
during the expansion for a BEC after N = 40 Bloch cycles. It takes more than 100 ms of expansion 
for the interference pattern to acquire full contrast. In general, we find that the contrast is improved 
when the horizontally confining beam Li is not switched off abruptly but is ramped down slowly 
within the first 55 ms of time-of-flight, reducing the horizontal expansion rate. However, this 
happens at the cost of some additional momentum broadening along the vertical direction. Our 
imaging techniques allow us to resolve structure in momentum space on a scale below O.lhk in a 
single shot absorption image. 
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To understand the interference structure and its evolution in time, we compute the total BEC 
wave function in quasi-momentum space for a phase evolution at the different lattice sites given 
by Equation (1), as detailed in the Methods section. Fig. 2b shows the interference pattern for our 
experimental parameters according to this simple model. The experimental results are qualitatively 
very well reproduced by the model when we reduce a mt by a factor of 0.9 compared to the value 
deduced from our experimental parameters. This scale factor accounts primarily for the fact that 
our simple model does not take into account any horizontal dynamics. In particular, switching off 
L 2 when starting the evolution leads to an excitation of a radial breathing mode in the horizon- 
tal plane, reducing the density at each site and modulating it in time. To a first approximation, 
rescaling of a- mt accounts for this. Nevertheless, the agreement between the experiment and the 
analytical model indicates that the dominant driving mechanism for the wave function spreading 
and interference is the nonlinear phase evolution. In particular, phase coherence is not lost, in 
contrast to previous experiments 22 . We test this coherence and demonstrate control over the phase 
evolution in two experiments. 

3 Cancellation of dephasing using an external potential 

Equation (1) suggests that the effect of interactions can be cancelled by the application of an ex- 
ternal potential. Indeed, choosing this potential to be equal to the initial loading potential, i.e. 
choosing a int « (3 tr , allows us to observe persistent Bloch oscillations for an interacting BEC. 
The BEC quasi-momentum distribution after N = AO Bloch cycles is shown in Fig. 4a and 4b as 
a function of the strength of the external compensating potential, given by the power in laser L 2 . 
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When the external potential does not compensate for interactions, the condensate wave function is 
dephased and spreads over the whole BZ within less than N = 20 Bloch cycles. In contrast, when 
the external potential balances the effect of interactions, the BEC wave function does not spread 
out and Bloch oscillations are clearly visible. The time during which Bloch oscillations can be ob- 
served is now greatly extended compared to the case when a compensating potential is absent. The 
transition from a dephased to a non-dephased wave function as a function of confinement strength 
is quantified in Fig. 4c, where the rms-width Ap of the singly-peaked quasi-momentum distribu- 
tion after N = AO Bloch cycles is plotted as a function of the laser power in L 2 . Fig. 4d and 4e 
show the time evolution of the quasi-momentum distribution without and with the compensating 
potential while all other parameters are kept the same. Fig. 4d essentially shows the broadening 
of the distribution as described before. Interestingly, the condensate wave function in the presence 
of a compensating potential shown in Fig. 4e dephases in a completely different way. Initially, the 
central peak shows no broadening. However, it is slowly depopulated, while a much broader back- 
ground distribution is increasingly populated. After about 100 oscillations, the shape of the central 
peak starts to develop side lobes or splits in two, with the exact shape varying from one experimen- 
tal run to the next. The timescale for the loss of interference is a factor 10 larger than the timescale 
on which the dephasing and hence the initial broadening takes place in the uncompensated case. 

4 Rephasing of a dephased condensate 

Second, we perform a matter wave spin-echo-type experiment. We initially proceed as shown in 
Fig. 2, letting the wave function evolve for a time corresponding to about N = AO Bloch cycles 
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until it is fully dephased and shows, upon measurement, a regular interference structure. We then 
essentially remove the effect of interactions by ramping to a = 10 a within 10 ms. By not switching 
the interaction entirely off and by ramping comparatively slowly we avoid excessive excitation of 
the radial breathing mode as a result of the change in the mean field potential at each site. At the 
same time, we gradually turn on the harmonic potential as given by the horizontal dipole trapping 
laser beam L 2 within 4 ms to approximately the same depth as during the initial BEC loading 
phase. From equation (2) we expect that the wave function now experiences a phase shift with 
a quadratic spatial dependence with opposite sign, allowing us to reverse the evolution and to 
recover the initial condition. Fig. 5 shows the resulting quasi-momentum distributions. As time 
progresses, the wave function indeed refocuses while it continues to perform Bloch oscillations. 
As we do not control the value of 8 for a particular run, we record about 10 distributions for each 
evolution time and select those that are symmetrical, corresponding to Bloch cycle phase = or 
4> = ir. For the chosen strength of the potential, refocusing happens after about 24 Bloch cycles 
after the ramp of a. This confirms that the initial broadening and dephasing mechanism must have 
been coherent. We note that we cannot avoid some excitation of the radial breathing mode as seen 
in the absorption images given in Fig. 5. 

5 Discussion 

Our results raise several important questions: To what extent can matter wave interferometry be 
performed in the presence of interactions? What sets the timescale for the eventual loss of inter- 
ference contrast? Certainly, our simple analytic model does not predict any loss of contrast. In 

10 



particular, it should be possible to completely eliminate the effect of interactions with the compen- 
sating external potential. However, there are several effects not included in the model that could 
cause the residual dephasing we observe. Motion in the radial direction, which causes the density 
and therefore the interaction energy to change over time, could lead to mixing of the different de- 
grees of freedom and hence to additional dephasing. This might apply to our matter wave spin-echo 
experiment, but in the experiment where we compensate by means of the external potential there 
is hardly any radial excitation and this effect should not play a role. The appearance of dynamical 
instabilities 23 25 can be ruled out, as our experimental parameters are outside the unstable region. 
Going beyond the mean-field treatment, a variety of factors can lead to dephasing. For example, at 
each lattice site there exists a superposition of number states, accumulating different phases corre- 
sponding to their respective interaction energies 26 - 27 . This leads to an effective dephasing, as the 
phase on a particular lattice site becomes ill-defined. Basic estimates 26 ' 27 indicate a dephasing time 
of about 130 ms for our system, on the same order as we observe. 

These experiments constitute a clear demonstration of coherent dynamics in an interacting 
macroscopic quantum system. This coherence affords a large degree of control over the system, as 
demonstrated by the possibility to rephase the wave function using an external potential in order to 
reverse dephasing due to interactions. The control demonstrated here has potential application in 
matter-wave interferometry, and such a degree of control over the mean-field evolution also opens 
the possibility to probe beyond-mean- field effects in atom interferometers. 



11 



Methods 



Sample preparation 

Our experimental approach initially follows the procedure described in ref. 12 . In brief, within 
10 s we produce an essentially pure BEC with tunable interactions 16 in the Thomas-Fermi limit 
with up to 1.5 x 10 5 Cs atoms. The BEC is trapped in a crossed-beam dipole trap generated by 
a vertically (Li) and a more tightly focused horizontally (L 2 ) propagating laser beam. The BEC 
is cigar-shaped with the long axis oriented along the direction of L 2 . The trap frequencies are 
(u x , uj y , uo z ) = 2jx (39, 5, 39) Hz, where x denotes the horizontal direction perpendicular to L 2 , 
y is the axial direction along L 2 , and z is the vertical direction. We magnetically control the 
scattering length a in the range between a and 300 a with a resolution of about 0. 1 a . For BEC 
production, we work at a = 210 a , where three-body losses are minimized 28 . Initially, we support 
the optical trapping by magnetic levitation against gravity 16 . As shown in Fig. la we superimpose 
an optical lattice with d = \/2 along the vertical direction, where A = 1064.5 nm is the wavelength 
of the lattice light. To load the BEC into the lattice, we stiffen the horizontal confinement within 
1 s, leading to trap frequencies of 2ir x (41, 13, 39) Hz, and at the same time turn on the lattice 
potential exponentially to a depth of 8E R . Here, E R = h 2 /(2m\ 2 ) = ks x 64 nK is the photon 
recoil energy and m the mass of the Cs atom. The BEC is thus gently loaded into the lattice, 
occupying about 25 to 35 lattice sites, with up to 7000 atoms at the central site. 
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Derivation of the BEC wave function in momentum space 

Here, we outline the method used to calculate the images in Fig. 2b. Due to the comparatively 
small interaction energies in our system, the atoms are restricted to move in the lowest Bloch band 
and we can write the local wavefunction at lattice site j as *ffj(r±, z) = Wo\z)$±(pj,r±), where 
Wq\z) is the lowest-band Wannier function localized at the j-th site and $±(rij, r±) is a radial 
wave function depending on the occupation number rij = \cj\ 2 at each site 17 . We can then write 
the total time-dependent wave function in momentum space as 

j 

= w^ ) ip z )Y J C j {t)e- ip ^i.{n j ,PA_). (2) 
j 

Transforming to quasi-momentum space and assuming that the phase at each lattice site evolves 
according to Equation (2), we can write 21 

V(q z ,p±,t) = ^c J (0)e~ i (' ?+ ^^e- l(/3trj2 - Qinu2) ' /?l $ ± (n j ,p ± ), (3) 

3 

where q z denotes the quasimomentum. The images in Fig. 2b show the BEC density distribution 
\^(Qz, P±, t) | 2 integrated along one radial direction, using a Thomas-Fermi wave function as radial 
wave function &±(rij,p ± ). 

We have compared the result in Fig. 2b with a numerical integration of the discrete nonlinear 
Schrdinger equation 17 , which includes tunneling between lattice sites, and find essentially identical 
results, confirming that tunneling is inhibited. 



13 



References 



1. Davisson, C. J., & Germer, L. H. The scattering of electrons by a single crystal of nickel. 
Nature 119, 558-560 (1927). 

2. Halban, H. v. Jr & Preiswerk, P. Preuve exprimentale de la diffraction des neutrons. C.R. Acad. 
Sci. 203, 73-75 (1936). 

3. Estermann, I. & Stern, O. Beugung von Molekularstrahlen. Z. Phys. 61, 95125 (1930). 

4. Andrews, M. R. et al. Observation of interference between two Bose condensates. Science 
275, 637-641 (1997). 

5. Berman, P. R. (ed.) Atom Interferometry (Academic, New York, 1997). 

6. Cronin, A. D., Schmiedmayer, J. & Pritchard, D. E. Atom interferometers. Preprint at 
< |http://arxiv.org/abs/07 12.3703> (2007). 

7. Hart, R. A., Xu X., Legere R. & Gibble, K. A quantum scattering interferometer. Nature 446, 
892-895 (2007). 

8. Wicht, A., Hensley, J. M., Sarajlic, E. & Chu, S. A preliminary measurement of the fine 
structure constant based on atom interferometry. Phys. Scr. T102, 82-88 (2002). 

9. Clade, P. et al. Determination of the fine structure constant based on Bloch oscillations of 
ultracold atoms in a vertical optical lattice. Phys. Rev. Lett. 96, 033001 (2006). 



14 



10. Fixler, J. B., Foster, G. T., McGuirk, J. M. & Kasevich, M. A. Atom interferometer measure- 
ment of the Newtonian constant of gravity. Science 315, 74-77 (2007). 

11. Gupta, S., Dieckmann, K., Hadzibabic, Z. & Pritchard, D. E. Contrast interferometry using 
Bose-Einstein condensates to measure h/m and a. Phys. Rev. Lett. 89, 140401 (2002). 

12. Gustavsson, M. et al. Control of interaction-induced dephasing of Bloch oscillations Phys. 
Rev. Lett. 100, 080404 (2008). 

13. Fattori M. et al. Atom interferometry with a weakly interacting Bose-Einstein condensate. 
Phys. Rev. Lett. 100, 080405 (2008). 

14. Ben Dahan, M., Peik, E., Reichel, J., Castin, Y. & Salomon, C. Bloch oscillations of atoms in 
an optical potential. Phys. Rev. Lett. 76, 4508^1511 (1996). 

15. Anderson, B. P. & Kasevich, M. A. Macroscopic quantum interference from atomic tunnel 
arrays. Science 282, 1686-1689 (1998). 

16. Weber, T., Herbig, J., Mark, M., Nagerl, H. & Grimm, R. Bose-Einstein condensation of 
cesium Science 299, 232-235 (2003). 

17. Smerzi, A. & Trombettoni, A. Nonlinear tight-binding approximation for Bose-Einstein con- 
densates in a lattice. Phys. Rev. A 68, 023613 (2003). 

18. Deng, L. et al. Temporal, matter-wave-dispersion Talbot effect. Phys. Rev. Lett. 83, 5407-541 1 
(1999). 



15 



19. Kaplan, A. E., Marzoli, I., Lamb, W. E. & Schleich, W. P. Multimode interference: Highly 
regular pattern formation in quantum wave-packet evolution. Phys. Rev. A 61, 032101 (2000). 

20. Kastberg, A., Phillips, W. D., Rolston, S. L., Spreeuw, R. J. C. & Jessen, P. S. Adiabatic 
cooling of Cesium to 700 nK in an optical lattice. Phys. Rev. Lett. 74, 1542-1545 (1995). 

21. Witthaut, D., Werder, M., Mossmann, S. & Korsch, H. J. Bloch oscillations of Bose-Einstein 
condensates: Breakdown and revival. 71, 036625 (2005). 

22. Morsch, O. et al. Decay and revival of phase coherence of a Bose-Einstein condensate in a 
one-dimensional lattice. Phys. Rev. A 67, 031603 (2003). 

23. Zheng, Y., Kostrun, M. & J. Javanainen. Low-acceleration instability of a Bose-Einstein Con- 
densate in an optical lattice. Phys. Rev. Lett. 93, 230401 (2004). 

24. Cristiani, M. et al. Instabilities of a Bose-Einstein condensate in a periodic potential: an ex- 
perimental investigation. Opt. Express 12, 4-10 (2004). 

25. Fallani, L. et al. Observation of dynamical instability for a Bose-Einstein Condensate in a 
moving ID optical lattice. Phys. Rev. Lett. 93, 140406 (2004). 

26. Li, W., Tuchman, A. K., Chien, H. & Kasevich, M. A. Extended coherence time with atom- 
number squeezed states. Phys. Rev. Lett. 98, 040402 (2007). 

27. Imamoglu, A., Lewenstein, M. & You, L. Inhibition of coherence in trapped Bose-Einstein 
condensates. Phys. Rev. Lett. 78, 2511-2514 (1997). 



16 



28. Kraemer, T. et al. Evidence for Efimov quantum states in an ultracold gas of caesium atoms. 
Nature, 440, 315-318 (2006). 

Acknowledgements We thank E. Arimondo, O. Morsch, W. Schleich, A. Smerzi, D. Witthaut and A. 
Buchleitner and his group for helpful discussions. We are indebted to R. Grimm for generous support and 
gratefully acknowledge funding by the Austrian Ministry of Science and Research (Bundesministerium fiir 
Wissenschaft und Forschung) and the Austrian Science Fund (Fonds zur Forderung der wissenschaftlichen 
Forschung) in form of a START prize grant and through SFB 15. R.H. is supported by a Marie Curie 
International Incoming Fellowship within the 7th European Community Framework Programme. 

Competing Interests The authors declare that they have no competing financial interests. 

Correspondence Correspondence and requests for materials should be addressed to H.-C. N. 
(email :christoph. naegerl @ uibk. ac. at) . 



17 



L, lattice 




Figure 1: BEC-based atom interferometer, a, Experimental configuration: The tunable BEC 
is formed at the intersection of the vertical guide laser beam Li and a horizontal trapping beam 
L 2 . The lattice is oriented along the vertical direction. Gravity, g, is initially compensated by a 
force due to a magnetic field gradient, VB. b, Imaging the first Brillouin zone (BZ): One cycle of 
Bloch oscillations for a non-interacting BEC as seen in time-of-flight absorption imaging, showing 
narrow peaks cycling through quasi-momentum space for cycle phases = 0, ix /A, it/2, to 2tt. 
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Figure 2: Interaction induced macroscopic matter wave interference, a, Experimental results 
showing the quasi-momentum distribution as a function of evolution time r given in units of the 
Bloch period. The absorption images are taken in steps of 4 Bloch cycles for a BEC with an initial 
peak density of n = 4 x 10 13 atoms/cm 3 loaded into about 35 lattice sites with a = 190 a . Each 
image corresponds to a single realization of the experiment, b, Evolution of the wave function in 
quasi-momentum space when the phase at the individual lattice sites evolves according to equation 
(2) with (3 tT = (no external trap) for n = 4 x 10 13 atoms/cm 3 loaded into 35 lattice sites with 
a = 190 a . «int is slightly rescaled to account for the reduction in density due to transversal 
dynamics, see text. In a, some additional broadening, largely due to the presence of the horizontal 
trapping potential during expansion, can be seen. 
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Figure 3: Contrast of interference fringes. Contrast of matter wave interference emerging during 
time- of- flight expansion for a BEC after N = 40 Bloch cycles, where the wave function completely 
fills the BZ. We define the contrast as (I max - Imin)/{I m ax + I m in), where I max (I min ) is the 
average value of the maxima (minima) of the central peak structure. Each data point is the average 
contrast of 10 experimental runs and the error bars indicate the la statistical error. The insets show 
measured quasi-momentum distributions integrated along the transverse direction at two expansion 
times as indicated. 
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Figure 4: Cancellation of interaction induced dephasing and observation of persistent Bloch 
oscillations, a-c, Absorption images showing the quasi-momentum distribution for cycle phase 
<p = tt (a) and = (b) after N = AO Bloch cycles and (c) momentum width Ap for = as a 
function of confinement strength, normalized to the confinement strength at loading, d Momentum 
distribution for <\> = as a function of the number N of Bloch cycles when no compensating 
potential is present, showing fast broadening, e The evolution of the momentum distribution for 
the case of optimum cancellation of interactions. 
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Figure 5: Matter wave spin-echo-type experiment. Rephasing of the BEC from a fully dephased 
wave function back into a narrow distribution after switching interactions to near zero and turning 
on an external potential. Time progresses from front to back. The black solid lines correspond to 
selected quasi-momentum distributions that refocus into the characteristic singly-peaked distribu- 
tion (cycle phase = 0), see text. They are separated in time by 1.15 ms or two Bloch cycles, and 
they are offset for clarity. The red solid lines correspond to selected distributions that refocus into 
the characteristic double-peaked distribution (cycle phase <p = tt). The images are absorption im- 
ages corresponding to the adjacent quasi-momentum distributions. Some radial excitation is also 
present. 
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